arXiv: 1502.05116vl [astro-ph.EP] 18 Feb 2015 


To APPEAR IN ApJ Letters 

Preprint typeset using style emulateapj v. 01/23/15 


THE FIRST SCIENCE RESULTS FROM SPHERE: DISPROVING THE PREDICTED BROWN DWARF 

AROUND V471 TAU* 

A. Hardy^’®, M.R. Schreiber^’®, S.G. Parsons C. Caceres^’®, G. Retamales^’® Z. Wahhaj^, D. MaweP, H. 

Ganovas^’®, L. Cieza®’®, T.R. Marsh®, M.C.P. Bours®, V.S. Dhillon,'^, A. Bayo®’® 

To appear in ApJ Letters 

ABSTRACT 

Variations of eclipse arrival times have recently been detected in several post common envelope binaries 
consisting of a white dwarf and a main sequence companion star. The generally favoured explanation 
for these timing variations is the gravitational pull of one or more circumbinary substellar objects 
periodically moving the centre of mass of the host binary. Using the new extreme-AO instrument 
SPHERE, we image the prototype eclipsing post-common envelope binary V47I Tan in search of 
the brown dwarf that is believed to be responsible for variations in its eclipse arrival times. We 
report that an unprecedented contrast of Amn = 12.1 at a separation of 260 mas was achieved, but 
resulted in a non-detection. This implies that there is no brown dwarf present in the system unless 
it is three magnitudes fainter than predicted by evolutionary track models, and provides damaging 
evidence against the circumbinary interpretation of eclipse timing variations. In the case of V471 Tau, 
a more consistent explanation is offered with the Applegate mechanism, in which these variations are 
prescribed to changes in the quadrupole moment within the main-sequence star. 

Keywords: stars: individual(V471 Tau) — planet-star interactions — binaries: eclipsing — brown 
dwarfs — white dwarfs — binaries: close 


1. INTRODUCTION 

Circumbinary substellar objects, although anticipated 
for a long time, have only rece ntly been identified around 
main-sequence binary stars (jPovle et al.l 1201111 . Long 
before this however, claims have been made for cir¬ 
cumbinary substellar objects around close white dwarf- 
main sequence binaries. Initially consisting of a main 
sequence binary with sepa rations of order ~1 AU 
(jZorotoyic fc Schreibeill201^ , these systems are belieyed 
to haye been through a spectacular phase of binary star 
eyolution to explain their current close separation. When 
the more massiye star of the binary eyolyes off the main 
sequence, it fills its Roche-lobe and generates dynami¬ 
cally unstable mass transfer onto the secondary star. As 
the time scale for dynamically unstable mass transfer 
is much shorter than the thermal time scale of the sec¬ 
ondary, the latter cannot adjust its structure fast enough 
to incorporate the oyerflowing mass. Instead, a common 
enyelope of material forms around both the secondary 
star and the core of the giant. Drag forces between 
the enyelope and the central binary then cause the enye¬ 
lope to be expelled at the expense of or bital energy and 
angular momentum of the binary (e.g. I Webbinfl 1 1 984 
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iZorotoyic et aLll^OlOHlyanoya et al.ll2013H . The emerging 
white dwarf-main sequence binaries contain separations 
of just a few solar radii, and are known as post common 
enyel ope binaries (PCEBs) (jNebot Comez-Moran et ^ 

l2onli . 

Shortly after the discoyery of the first PCEB it was 
realised that it displays yariations in its eclipse ar- 
riyal times. Today, similar yariations are seen in al- 
most all eclipsing PCEBs with long enough cov erage 
(IParsons et al.l I2010IJ : IZorotovic fc Schreibeii I2013I1 , for 
which the most common hypothesis is the presence of 
a circumbinary object, typically a brown dwarf or mul¬ 
tiple giant planets. In this scenario, the gravitational 
pull of the circumbinary objects periodically move the 
center of mass of the host binary stars, thereby chang- 
ing the ligh t travel time of the eclipse signal to Earth 
(|IrwinllT9^ . Indeed, the planetary model employed to 
explain the ecli pse timing variations (E TVs) seen in the 
PCEB NN Ser (|Beuermann et ani2010ll successfully pre¬ 
dicted new eclipse arrival times ( Beuermann et ani2013t 
iMarsh et al.l 1201^ . providing support to the circumbi¬ 
nary interpretation but raising q uestions regarding the 
forma tion of these third objects. IZorotovic fc Schreiberl 
(|20I3f) favour a scenario in which the circumbinary ob¬ 
jects form as a consequence of the common envelope evo¬ 
lution, in a so-called second generation scenario. This is 
based on the finding that nearly all PCEBs with long 
enough coverage show ETVs, yet only a small fraction of 
main sequence binarie s seem to host circumbinary sub¬ 
stellar objects. Indeed. ISchleicher fc DreizI^ (|2014[ 1 were 
able to develop a model in which a second generation 
protoplanetary disk forms during common envelope evo¬ 
lution and produces giant p lanets through the d isk in¬ 
stability model. In contrast. iBear fc Sokeil (I2014D prefer 
the first generation scenario, in which the objects form 
at a similar time to their main-sequence hosts, and sur¬ 
vive the common-envelope phase. They claim that if a 
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second generation scenario were true, too large a fraction 
of the common envelope mass would have to form into 
substellar companions. 

However, before investigating possible formation sce¬ 
narios further, we must exercise caution with the third 
body hypothesis. Although the circumbinary object 
model has proved successful in the case of NN Ser, this 
is an exception. In general the predictions from pro¬ 
posed planetary systems around PCEBs disagree with 
more recent eclipse t iming measurements (jParsons et al.l 
l2010bl: iBours et al.l 1201'^ . and some of the proposed 
planetary sy stems are dynamically unstable on very short 
time scales (|Wittenmver et al]l2013HHorner et'al]l2012L 
1201311 . The failure of all circumbinary object models ex¬ 
cept the one for NN Ser implies either that our timing 
coverage is insufficient, or that there must be an alterna¬ 
tive mechanism driving ETVs. 

To progress with this situation, it has become vi¬ 
tal that the circumbinary companion interpretation be 
tested independently. The most conclusive way to 
achieve this is to image one of the proposed objects and 
the natural choice for such an observation is V471 Tau. 
V471Tau consists of a 0.84T0.05 Mfn white dwarf an d 
a O.93±O.O7M0 secondary star (|0 IBrien et al.l 120011 1. 
and is a member of the 62 5 Myr old Hyades open 
clust er (fPerrvman et al.l 1199811. Soon after its discov¬ 
ery (|Nelson fc Yound 1197011 . iLohsenI (1197411 reported 
ETVs which have been interpret ed as being caused 
by a circumbinary br own dwarf (iBeavers et al.l 119861 
Icininan fc RibasI 1200111 . V471Tau is ideal to test the 
circumbinary interpretation because it is nearby, bright, 
and the proposed brown dwarf reaches projected sep¬ 
arations exceeding 200 mas, making detection possible 
with the new extr eme-AO facilities such as SPHERE 
(jBeiizit et al.ll200'^ . 

Here we present new high-precision eclipse times of 
V471 Tau, and use these to refine the proposed brown 
dwarf parameters using the Markov Chain Monte Carlo 
(MCMC) method. We then test the circumbinary in¬ 
terpretation of ETVs with SPHERE science verification 
observations, with high enough contrast to detect the 
brown dwarf independent of if it formed in a second or 
first generation scenario. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. High-speed eclipse photometry 

In order to refine the orbit of the circumbinary com¬ 
panion we obtained high-speed photometry of the eclipse 
of the white dwarf in V471 Tau with the frame-transfer 
camera ULTRACAM (iDhillon et al.l[200^ mounted as a 
visitor instrument on the 3.6-m New Technology Tele¬ 
scope (NTT) on La Silla in November and December 
2010. ULTRACAM uses a triple beam setup allowing 
one to obtain data in the u', g' and i' bands simultane¬ 
ously, with negligible dead time between frames. Due to 
the brightness of the target we de-focussed the telescope 
and used exposures times of ^5 s. There were no suit¬ 
ably bright nearby stars in ULTRACAM’s field of view 
to use for comparison. We therefore used the i' band 
observations, where the eclipse is extremely shallow, as 
a comparison source for the u' band data. This results 
in a large colour term in the resulting light curve, but 
does remove variations in the conditions and does not 


Table 1 

ULTRACAM eclipse times for V471 Tau 


Cycle 

HMJD(TDB) 

Uncertainty (seconds) 

2848 

5512.2840584 

1.76 

2886 

5532.0889885 

1.59 

2911 

5545.1185942 

1.62 

2915 

5547.2033608 

2.37 


affect the eclipse timings. In some observations the con¬ 
ditions were good enough that no comparison source was 
required. 

These data wer e reduced us i ng th e ULTRACAM 
pipeline software (|Dhillon et al.l 120071) and the resul¬ 
tant eclipse light curves were fitted with a dedicated 
code designed to fit bina ries containing white dwarfs 
(|CoDDerwheat et ^I2010H . The measured eclipse times 
were then heliocentrically corrected and are listed in Ta¬ 
ble. 1. 

2.2. SPHERE observations 

The imaging data testing the existence of the brown 
dwarf were acquired using the extreme adaptive optics 
instrument, SPHERE at the UTS Nasmyth focus of the 
VLT, on 2014 December 11. An earlier set of obser¬ 
vations was performed on 2014 December 6, but is not 
used here because of poorer data quality. SPHERE is 
able to provide H and K-band images with Strehl ratios 
> 90%. The integral field spectrograph (IFS) and in¬ 
frared dual-band imager and spectrograph (IRDIS) were 
used simultaneously in the IRDIFS mode. The IFS de¬ 
livered a datacube containing 38 monochromatic images 
at a spectral resolution of R'^50 between 0.96-1.34 ^m, 
whilst IRDIS delivered dual-band imaging in the H2 and 
H3 bands (central wavelengths of 1.59 /im and 1.67/rm 
respectively, and FWHMs of 0.0531/J,m and 0.0556/rm). 

The brightness of the target and desired contrast re¬ 
quired the use of the N_ALC_YJH coronagraph, with 
an inner working angle of 0.15”. Detector integration 
times were set at 64 s for each detector. The observa¬ 
tions were obtained in pupil-stabilized mode, where the 
held rotates. To derotate and combine the images, one 
needs to accurately measure the center of rotation which 
is also the location of the star behind the coronagraph. 
To achieve this, a waffle pattern was introduced into the 
deformable mirror of the AO system, placing four repli¬ 
cas of the star equidistant from it in a square pattern. 
These calibration images were taken before and after the 
science, allowing subpixel accuracy in centroiding. Off- 
coronagraph, unsaturated images of the star were also 
obtained with a neutral density hlter to allow sensitiv¬ 
ity/contrast measurements. Basic reduction, including 
dark and hat-helding was performed using the SPHERE 
pipeline. Custom angular differential imaging (ADI) rou¬ 
tines, particularly for subpixel shifting a nd aligning of 
speck les, were used to subtract the star (iWahhai et al.l 
l2013f) . A custom principal component analysis routine 
was also a pplied to the data to compare with the ADI 
reduction (IMawet et al.ll201§l . 

3. PREDICTING THE PROJECTED SEPARATION OF THE 
POTENTIAL BROWN DWARF 
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Assuming a third body orbiting around V471 Tau, the 
time delay or advance caused by this body can be ex¬ 
pressed as 


AT = 


012 sin is 
c 


1 - e§ 


1 - 1-63 cos P 3 


sin(p 3 + bj) + es sin(tLi) 


(e.g. lirwini [195^ where c is the speed of light and 032 
is the semi-major axis of the binary star’s orbit around 
the common center of mass of the triple system. The 
other parameters define the orbit of the third body, i.e. 
its inclination is, the orbital eccentricity 63 , argument of 
periastron w, a nd true anom a ly vs- 

As shown by IMarsh et al.l (120141) , strong correlations 
can exist between orbital parameters and the problem 
is highly degenerate unless a large number of high pre¬ 
cision eclipse timing measurements are available. Only 
our recent ULTRACAM measurements provide precise 
eclipse timings, with uncertainties of ~ 1 . 8 s, whereas the 
timings in the literature have been assigned large error 
estimates of 15 s for the sake of caution. To properly 
identify not only the best fit parameters but also their 
uncertainties we performed a Markov Chain Monte-Carlo 
(MCMC) simulation to the eclipse times. 

The prediction of the best-fit model can be seen in 
Fig. 1, (top left panel) alongside all archival observed- 
minus-calculated eclipse times (|Kundra fc Hrid 120111) 
and the new times reported in this paper. This best-fit 
model corresponds to a brown dwarf of mass 0.044±0.001 
Mq and semi-major axis 12.8±0.16 AU. 

While the parameters of the brown dwarf in the one 
body fit are well constrained by the eclipse times (Fig. 2), 
the residuals are far from random and suggest another 
mechanism may also be at work. To test this possibil¬ 
ity, we performed another MCMC with 2 companions to 
account for these deviations. This further allowed us to 
test whether the brown dwarf causing the main variation 
could be at a smaller separation or be less massive, which 
would make it harder to detect. The resulting best fit is 
shown in the bottom left panel of Fig.[TJ 

The derived orbital parameters for both cases were 
then projected onto the sky using a distanc e to V471 Tau 
of_~50 pc, as measured by Hipparcos (Ivan LeeuwenI 
I2007D . to obtain a predicted separation for the brown 
dwarf companion in December 2014. The simulation sug¬ 
gests a separation of 260f(fgmas for the one body ht and 
this value hardly changes if an additional object is as¬ 
sumed to account for the problems of the one body ht 
(see Fig. 1, right panels). 


4. TESTING THE PREDICTION WITH SPHERE 

No third component is present in the SPHERE IRDIS 
images (Fig. 3, left panel). The contrast achieved was 
estimated via two different methods of fake companions 
injection. In the hrst, fake companions of a known con¬ 
trast were injected at different angular separations, and 
the contrast dehned by wh ere the fake compan ion was 
recovered 95% of the time (|Wahhai et al.ll2?)l^ . In the 
second method, the fake companions were used to mea¬ 
sure the post-ADI throughput loss, and this was used 
to renormalize the contrast curve of a typical saturated 
SPHERE PSF. In both cases, the H2 and H3 channels 
were summed as no spectral difference between the chan¬ 
nels was expected, and the curves were corrected for 


small -sample statistics at small separations (jMawet et al.l 
120141) . The resulting contrast curves for the IRDIS de¬ 
tector can be seen in Fig. 3, right panel, with the former 
method as the solid line, and the latter as the dashed. 
There is good agreement between the 2 methods at the 
predicted separation of ^260 mas, and both indicate an 
achieved contrast of ~12.1 magnitudes. To determine 
if this is sufficient to indeed detect the brown dwarf, it 
is necessary to know its age, mass and metallicity, from 
which the brown dwarf luminosity can be predicted. We 
find the mass is well constrained from the MCMC mod¬ 
els, and the metallicity was assumed identical to other 
members of the Hyade s cluster, with [M/H]=0.14±0.05 
(jPerrvman et al.lH998D . 

The cooling age of the white dwarf in V471 Tau is 
~ lOMyrs, and puts a stringent constraint on the age 
in a second generation scenario. If the 0.044 Mq brown 
dwa rf had formed in such a scenario, BTSettl mod- 
els (lAlla rd et al l 1201211 combined with isocrones from 
iBaraffe et al.l (I1998D predict a Amu ~ 4.5. This is 7.5 
magnitudes brighter than our detection limit, conclu¬ 
sively ruling out a second generation formation scenario 
for the potential brown dwarf around V471 Tau. 

If the brown dwarf formed in a first generation sce¬ 
nario, we can obtain an estimate of its age from the age 
of the Hyades cluster (625 Myr). An identical modelling 
procedure suggests that such a brown dwarf will have 
a contrast of Amn ~ 9.2 in the H band, 3 magnitudes 
brighter than our detection limit. 

5. DISCUSSION 

5.1. Could the BD have escaped our SPHERE 
observations? 

Inspecting the contrast curve presented in Fig. 3, it 
is clear that if a first generation brown dwarf was at a 
smaller separation and/or had a lower mass, it may still 
have escaped detection. To test whether a brown dwarf 
could escape detection whilst simultaneously causing the 
0-C variations, we performed a final MCMC simulation 
with separations limited to only those which would be 
undetectable by SPHERE. The resulting fit drastically 
failed to explain the 0-C data, and delivered a of 
47.7, compared to 1.1 in the case of the constraint-free 
brown dwarf. As shown in Fig. 1, (lower left panel), not 
even the addition of an extra body can cause the brown 
dwarf to reach a separation undetectable by SPHERE. 
In order to recreate the observed 0-C variations, a large 
mass, high separation companion seems needed, and we 
therefore find no configuration in which the brown dwarf 
would have escaped SPHERE. 

The second possibility as to how the brown dwarf 
might have escaped detection, is that the models of sub- 
stellar evolutionary tracks are incorrect. Indeed, if the 
0.044 M 0 brown dwarf was at a temperature of ^700 K, 
as oppose to the ^1500 K predicted for a hrst-generation 
brown dwarf, it would escape detection. However, de¬ 
spite the uncertainties that are associated with these 
cooling models, th ere is no evidence to suggest models a re 
off to this degree (jBell et al.ll2012l : ILodieu et al.ll2014D . 

Finally, the brown dwarf around V471 Tau might be 
unique because of the evolutionary history of the host 
binary star. One might for example speculate that it ac¬ 
creted significant amounts of material during the rather 
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Figure 1. Observed-minus-calculated (0-C) eclipse arrival times for V471 Tau. Upper left panel: The eclipse times and associated errors, 
with the high precision eclipse times reported in this paper as the red open circles. The best-fit model assuming a third-body perturbation 
shown as the solid black line. The residuals of this fit can be seen directly beneath the curve. Upper right panel: The projected separations 
predicted by the MCMC simulation for observation in December 2014, assuming a distance of 50 pc to V471 Tau. The lower panels are 
identical, but denote the results of a 2 body fit the the eclipse timing variations. The dashed lines in the lower left panel denote the 
contributions from the different objects, with their sum in black. Despite the extra body, the predicted separation and mass of the brown 
dwarf hardly change. 



Figure 2. Correlations and histograms for the Markov-chain Monte Carlo simulation. M is the mass of the brown dwarf, a the semi-major 
axis of its orbit, Pq the period, e the eccentricity, oj the argument of periastron, the time of periastron of the binary, Psin the period 

of the binary, and Tq the time of periastron of the brown dwarf. All parameters appear well constrained with no significant correlations. 
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RA(mas) 



Separation(mas) 


Figure 3. H-band image of V471 Tau obtained on the SPHERE IRDIS instrument at the VLT. Left panel: Resulting image after angular 
differential imaging (ADI). The area in-between the white circles denotes the 5 sigma predicted position of the brown dwarf. Right panel: 
Contrast curves obtained via 2 different methods of fake companion injection (see text for details). The vertical lines again denote the 
predicted position of the object, with the diamond denoting the predicted contrast of a first generation brown dwarf. 
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recent common envelope phase, and this could have 
caused its metallicity to deviate significantly from the 
metallicity of the Hyades. Indeed, if the brown dwarf 
were hrst generation, and possessed a metallicity [M/H] 
= 0.5, models predict it would not be det ected in the pre¬ 
sented observations (|Allard et al.ll2niTl i. However, the 
metallicity of the secondary K star has been studied, 
and found to show no peculiar abundances despite the 
possibility i t accreted material du ring the common enve¬ 
lope phase (jStill fc Hussainll^OO.Hi , It is therefore highly 
unlikely that the brown dwarf metallicity is high enough 
to explain the non-detection. A final effect of recent ac¬ 
cretion, could be that the infallin g mass made the br own 
dwarf fainter due to contraction ( Baraffe et al.l[200^ . al¬ 
though confirming if this is indeed possible will require 
detailed evolutionary brown dwarf models dedicated to 
V471 Tau which is beyond the scope of the present paper. 

5.2. Applegate’s mechanism 

If our non-detection is due to the lack of a brown dwarf, 
an alternative mechanism must be responsible for the 
ETVs. Currently, there are two such alternative theories 
for ETVs in PCEBS. The first, apsidal precession, can 
confidently be ruled out for V471 Tau. This mechanism 
prescribes ETVs to the time dependance in the argu¬ 
ment of periastron. Although apsidal precession would 
require an eccentricity in the binary of just e=0.01 to 
create the ~300s timing variations seen in Fig.[l] the 
timescale would be less than 3 years, and not ^30 years 
as observed. The sec ond alternative e xplanation, the Ap¬ 
plegate mechanism ()ADDleefatelll992[l . is potentially able 
to drive the variations of the eclipse arrival times seen 
in V471 Tau. This theory prescribes these variations to 
quasi-periodic oblateness changes in the main sequence 
star, presumed to be driven by solar-like magnetic cy¬ 
cles. The K star in V471 Tau is particularly active, and 
may provide sufficient energy to drive these variations. 
Applegate’s mechanism is therefore the currently most 
convincing explanation for the eclipse arrival times ob¬ 
served in V471 Tau. 

However, although a suitable explanation in the case 
of V471Tau, in almost all other PCEBs showing ETVs, 
classical Applegate’s mechanism can be ruled out as they 


tend to contain less active main-sequence stars. NN 
Ser is one such system, and currently onl y the proposed 
planetary system can explain i ts behavior (iParson s et al.l 
12014 iBrinkworth et al.l 120061 : IParsons et al.ll2010all . al- 
though it is possible that an as-yet unconsidered model 
of magnetic field variations could still offer an explana¬ 
tion. Thus, Applegate does not offer a comprehensive 
explanation for ETVs, and confirmation of its effect in 
V471 Tau is needed alongside additional tests of the third 
body interpretation in other systems. 

6. CONCLUSION 

We have presented deep SPHERE science verification 
observations of V471 Tau testing the hypothesis that the 
observed ETVs are caused by a circumbinary brown 
dwarf. We reached an excellent contrast of Amn = 12.1 
at the predicted separation of the brown dwarf but no 
companion can be seen in the images. This excludes both 
a brown dwarf formed in a second generation scenario, as 
well as a standard brown dwarf at the age of the Hyades 
cluster to be present around V471 Tau. The Applegate 
mechanism is hence the one and only remaining model 
currently explaining the ETVs seen in V471 Tau. 

With this result, the origin of ETVs in PCEBs remains 
a puzzle. While no theory but the existence of two cir¬ 
cumbinary planets can currently explain the variations 
seen in the PCEB NN Ser, the most reasonable expla¬ 
nation for the variations seen in V471Tau is now the 
Applegate mechanism. We therefore conclude that in 
their current form neither the third body interpetation 
nor Applegate’s mechanism offer a general explanation 
for the ETVs observed in nearly all PCEBs. 
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